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Abstract. We present the new results of the Wuppertal-Budapest lattice QCD 
collaboration on flavor diagonal and non-diagonal quark number susceptibilities with 
2+1 staggered quark flavors, in a temperature regime between 120 and 400 MeV. 
A Symanzik improved gauge and a stout-link improved staggered fermion action is 
utilized; the light and strange quark masses are set to their physical values. Lattices 
with Nt = 6, 8, 10, 12 are used. We perform a continuum extrapolation of those 
observables for which the scaling regime is reached, and discretization errors are under 
control. 
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Figure 1. Left: Diagonal light quark susceptibility as a function of the temperature. 
Right: Non-diagonal u — s susceptibility as a function of the temperature. In both 
panels, data are obtained on Nt — 6, 8, 10 and 12 lattices. Two different spatial 
volumes are considered for the latter: 32'^ and 36'^. The band in the left panel is the 
comparison with the Hard Thermal Loop predictions taken from Ref. |11| 

1. Introduction 

A transition occurs in strongly interacting matter from a hadronic, confined system at 
small temperatures and densities to a phase dominated by colored degrees of freedom 
at large temperatures or densities. This field of physics is particularly appealing 
because the deconfined phase of QCD can be produced in the laboratory, in the 
ultrarelativistic heavy ion collision experiments at CERN SPS, RHIC at Brookhaven 
National Laboratory, ALICE at the LHC and the future FAIR at the GSI. A lot of effort 
is invested, both theoretically and experimentally, in order to find observables which 
would unambiguously signal the QCD phase transition. Correlations and fluctuations 
of conserved charges have been proposed long ago to this purpose [H |2]. The idea is 
that these quantum numbers have a very different value in a confined and deconfined 
system, and measuring them in the laboratory would allow to distinguish between the 
two phases. 

Fluctuations of conserved charges can be obtained as linear combinations of 
diagonal and non-diagonal quark number susceptibilities, which can be calculated on the 
lattice at zero chemical potential [HI H]. In the present contribution we show the results 
of our collaboration on some of these observables, with 2+1 staggered quark flavors, in 
a temperature regime between 120 and 400 MeV. The light and strange quark masses 
are set to their physical values. Lattices with A^^ = 6, 8, 10, 12 are used. Continuum 
extrapolations are performed for those observables for which discretization errors are 
under control and the scaling regime is reached. 

2. Results 

Quark number susceptibilities are defined as derivatives of the QCD pressure with 
respect to the quark chemical potentials. In particular, we consider the following ones: 
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Figure 2. Left: Diagonal strange quark susceptibility as a function of the 
temperature. Data are obtained on Nt = 6, 8, 10 and 12 lattices. Two different 
spatial volumes are considered for the latter: 32'^ and 36^. The red band indicates the 
continuum extrapolation. The black curve is the Hadron Resonance Gas (HRG) model 
prediction. Right: comparison between the light and strange quark susceptibility as 
functions of the temperature, for Nt = 12. 

The lattice action is the same as we used in [5l [6] , namely a tree-level Symanzik improved 
gauge, and a stout-improved staggered fermionic action (see Ref. [7] for details). The 
stout-smearing p] reduces the taste violation: this kind of smearing has one of the 
smallest taste violation among the ones used so far in the literature, together with the 
hisq action recently proposed by the hotQCD collaboration [9l [10]. In the left panel of 
Fig. [T] we show the diagonal light quark susceptibility as a function of the temperature. 
Results are obtained on five different lattices: Nt = 6, 8, 10 and two different spatial 
volumes for A^^ = 12: 32'^ and 36^. This observable reaches approximately 90% of the 
Stefan-Boltzmann limit already at T ~ 400 MeV. Notice that we do not perform a 
continuum extrapolation: this observable is pion-dominated at small temperatures and 
we need either finer lattices or an action which reduces even more the taste violation, 
in order to reach the scaling regime and give the correct continuum estimate. In the 
right panel of Fig. [1] we show the non-diagonal u — s susceptibility, which measures the 
correlation between these different flavors. As we can see, this observable has a peak in 
the vicinity of the phase transition, and tends to zero for infinitely large temperatures. 
This observable is very noisy and before performing a continuum extrapolation we need 
to reduce the errorbars for the Nt = 12 sets of data. In the left panel of Fig. |2] we 
show the strange quark number susceptibility as a function of the temperature. For this 
observable we perform a continuum extrapolation, as well as a comparison to the HRG 
model predictions. In the right panel of Fig. [2] we show, for A''^ = 12, a comparison 
between light and strange susceptibilities. We notice that the rapid rise as a function of 
the temperature occurs at larger temperatures for strange quarks than for light quarks. 

Baryon number susceptibility is defined as a linear combination of quark number 
susceptibilities in the following way: 

XB = ^ (2cr + + 2cf + 4cr' ■ ' TdHnZ 
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Figure 3. Left: Baryon number susceptibility as a function of the temperature. 
Right: Baryon-strangeness correlator as a function of the temperature. In both panels, 
results for Nt = 6, 8, 10 and 12 lattices are shown. Two different spatial volumes are 
considered for the latter: 32'^ and 36'^. The red band is the continuum extrapolated 
result, the black curve is the prediction of the HRG model. 

We show our results for this observable in the left panel of Fig. [3l a continuum 
extrapolation is performed, and compared to the HRG model results. The agreement is 
very good. In the right panel of Fig. |3] we show the baryon-strangeness correlator that 
was proposed in Ref. [12] as a diagnostic of strongly interacting matter. It is defined 
as: 
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The continuum extrapolation is in good agreement with the HRG model prediction for 
temperatures smaller than T^. 
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